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Abstract 

The process e"'"e~ W~'"W~7 is analysed using the data collected with the 
L3 detector at LEP at a centre-of-mass energy of 188.6 GeV, corresponding to an 
integrated luminosity of 176.8 pb~^. Based on a sample of 42 selected W"'"W~ 
candidates containing an isolated hard photon, the W"'"W~7 cross section, defined 
within phase-space cuts, is measured to be: crww7 = 290 ± 80 ± 16 fb, consistent 
with the Standard Model expectation. 

Including the process e"'"e~ — > z/z/77, limits are derived on anomalous contribu- 
tions to the Standard Model quartic vertices W"'"W~77 and W"^W~Z7 at 95% CL: 

-0.043 GeV-2 < ao/A^ < 0.043 GeV'^ 

- 0.08 GeY-^ < aJA^ < 0.13 GeV-^ 

- 0.41 GeV"^ < a„/A^ < 0.37 GeV'^ 



Submitted to Phys. Lett. B 



Introduction 



The LEP centre-of-mass energy for e'^e" collisions is now well above the kinematic threshold for 
W-pair production allowing for the study of radiative W-pair production, e"'"e~ W"'"W~7. 
The Standard Model (SM) [T1I2] predicts the existence of quartic gauge couplings (QGCs), 
leading to ^^^^^^^7 production via s-channel exchange of a 7 or Z boson as shown in Figure Ct- 

As the contribution of these two quartic Feynman diagrams with respect to the other com- 
peting diagrams, mainly initial-state radiation, is negligible at the LEP centre-of-mass energies, 
the process leading to the W"^W~7 final state could thus be sensitive to anomalous contribu- 
tions to the SM quartic gauge-boson vertices W+W~77 and W~^W~Z'y. 

The existence of Anomalous QGCs (AQGCs) would also affect the e"'"e~ i^ei^ell process 
via the W"'"W~ fusion Feynman diagram containing the W"'"W~77 vertex ^ (see Figure^). 
In the SM the reaction e+e~ uu'yy proceeds predominantly through s-channel Z exchange 
and t-channel W exchange, with the two photons coming from initial state radiation, whereas 
the SM contribution from the W^W~ fusion is negligible at LEP. AQGCs would enhance the 
uu'j'j production rate, especially for the hard tail of the photon energy distribution and for 
photons produced at large angles with respect to the beam direction. 

Here we describe the cross section measurement for the process e"'"e~ — *• W"'"W~7 and the 
determination of AQGCs using the data collected in 1998 with the L3 detector ^ at ^/s = 
188.6 GeV (denoted as y/s = 189 GeV hereafter) corresponding to an integrated luminosity 
of 176.8 pb~^. AQGCs are also independently determined using acoplanar photon-pair events 
with missing energy. This analysis is performed using the data at ^/s = 189 GeV and at 
y/s = 182.7 GeV collected in 1997 (denoted as y/s = 183 GeV hereafter) corresponding to a 
total integrated luminosity of 231.7pb~^. The results derived on AQGCs from the W"'"W~7 
and z/^'77 channels are finally combined. 

The search for anomalous contributions to the SM quartic couplings is performed within the 
theoretical framework of References El and El Recently, experimental measurements for these 
couplings have already been performed on final states with three vector bosons W"''W~7 [Z] 
and Z77 j8j. 



W+W 7 Final State and Signal Definition 

There are 14 Feynman diagrams at the tree level leading to the W^W~7 final state, and many 
other diagrams corresponding to photons from the decay products of hadronic or leptonic W's. 
We are interested only in two of these, the quartic diagrams. The other diagrams leading to 
the same final state are initial state radiation (ISR), final state radiation (FSR), and radiation 
from the W boson itself. 

The Monte Carlo used for the W"'"W~7 cross section determination is KORALW p. This 
generator does not include the quartic coupling diagrams. Initial state multi-photon radiation 
is implemented in KORALW in the full photon phase space. FSR from charged leptons in 
the event up to double bremsstrahlung is included using the PHOTOS [lOj package. Frag- 
mentation processes of quarks into hadrons are made according to the JETSET [JJj algorithm 
including photons in the parton shower. For the W~'"W~7 cross section measurement, this 
modelling is sufficient since the contribution of all the other diagrams is very small. The back- 
ground processes such as e+e" — > Z/7 — > qq(7) and e"'"e~ ZZ —> 4/(7) are simulated using 
PYTHIA [12\. The L3 detector response is simulated by the program GEANT jl3j . 
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In this analysis, the W"^W~7 signal is defined by the following phase-space cuts: 

• > 5 GeV, where is the energy of the photon, 

• 9^ > 20°, where 9^ is the angle between the photon and the beam axis, 

• > 20°, where is the angle between the direction of the photon and that of the 
closest charged lepton or jet. 

These cuts are mainly chosen for experimental reasons, to optimise the photon identification 
and the background suppression. They also largely avoid any infrared and coUinear singularities 
in the calculation of the signal cross section. 

The theoretically predicted W"'"W~7 cross section from KORALW corresponds to 272 ± 
4(stat) fb. In the EEWWG program jH] the effect of undetected additional ISR collinear to the 
beam pipe is included ^3] by implementing the EXCALIBUR ^H] collinear radiator function. 
The effect of the higher order radiative corrections is to move the effective centre-of-mass energy 
towards lower values, reducing the expected signal cross section by about 18%. The resulting 
EEWWG cross section corresponds to 233 ± 12(theor) fb. This is used as the SM expectation 
in the anomalous coupling analysis, which leads to less stringent constraints on AQGCs. The 
theoretical uncertainty ^3] is propagated to the AQGC determination. Consistent results are 
obtained with the YFSWW3 [IE] MC which predicts 224 ± 6(stat) fb. Differences of this 
order in the predicted cross section with high transverse momentum photons are expected J7] 
between pure leading-log and leading-log plus matrix-element based calculations. 

W+W~7 Event Selection and Cross Section 

The W-pair event selections used here are similar to those reported in Reference ^1 Only the 
semileptonic and fully hadronic W-pair decay modes are considered. The number of selected 
data events and the expected number of signal and background events are shown in Table 



Decay Channel 


Nobs 




i\tot 




qqei/e 


355 


0.768±0.005 


361 


19 




364 


0.834±0.002 


375 


19 




313 


0.605±0.003 


300 


42 


qqqq 


1514 


0.892±0.006 


1486 


296 



Table 1: Number of observed events, selection efficiencies with statistical uncertainties, ex- 
pected total number of events and background estimates for the various W"*'W~ decay channels 
according to the SM prediction. The efficiencies shown here include the contribution of cross 
efficiencies from the other W-pair decay modes. 

The photon selection in W"*"W~ events is optimised for each four-fermion final state. Photons 
are selected by requiring energy deposition in the electromagnetic calorimeter not associated 
with any track in the central detector, and low hadronic activity in a cone of half-opening angle 
of 7° around the electromagnetic cluster. The profile of the shower must be consistent with that 
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of an electromagnetic particle. In addition, the highest energy photon has to satisfy the signal 
definition requirements: > 5 GeV, 6*^ > 20°, and a-^ > 20°. Figure El shows the distributions 
of E^, 6^, and a-y, where for the last variable the direction of the reconstructed hadronic jet is 
assumed as the direction of the quark in the final state. In general, good agreement between 
the data and the SM expectation is observed. 

Table 121 summarises the results for the applied selection criteria. In total 42 W'''W~7 events 
are selected, where 37.8±0.6 is the Monte Carlo expectation. The efficiency £:ww7 is defined as 
the number of selected KORALW events (regardless of any phase-space cuts) divided by the 
number of generated MC events satisfying the signal definition. It accounts for small possible 
migration effects of events from outside the signal region into the selected sample due to the 
finite detector resolution. 



Decay Channel 


Nobs 






nrexp I Arexp 
^^FSR "T ^^Bkgr 


qqei/e7 


6 


0.483±0.025 


5.85±0.26 


2.26±0.14 




5 


0.547±0.027 


6.87±0.28 


2.88±0.17 


qqrz/^7 


7 


0.351±0.018 


4.63±0.22 


2.05±0.14 


qqqq7 


24 


0.504±0.016 


20.4±0.38 


9.23±0.26 


Total 


42 




37.8±0.6 


16.4±0.4 



Table 2: Number of observed events, selection efficiencies with statistical uncertainties and 
expected number of total and background events including final state radiation. 

The W"'"W~7 cross section is evaluated channel by channel and then combined according 
to the SM W-pair branching fractions. The result is: 

aww7 = 290 ± 80 ± 16 fb , 

where the first error is statistical and the second systematic. The measurement is in good 
agreement with both the KORALW and EEWWG SM expectations. 

Figure El shows the result obtained together with the predicted total W"''W~7 cross section 
from the EEWWG Monte Carlo as a function of the centre-of-mass energy. 

The systematic uncertainties arising in the inclusive W-pair event selections ^Hj are prop- 
agated to the final measurement and correspond to an uncertainty of ± 6.3 fb. Other possible 
systematic biases due to detector effects such as electro-magnetic cluster resolution, angular 
resolution, and calorimetric energy scale uncertainty are found to have a negligible effect on 
the final result. 

The total systematic uncertainty is dominated by the JETSET modelling of photons from 
meson decays (vr",?]). To estimate this effect, a data sample of 3.9 pb~^ collected in 1998 at 
^/s = 91 GeV is studied. The same photon identification criteria are applied to the selected 
Z — > qq events. An overall excess of (20 ± 10)% in the photon rate is found in data with 
respect to the PYTHIA Monte Carlo which uses the same JETSET fragmentation algorithm 
and particle decays as KORALW. A correction factor given by the ratio of photon production 
rates in data and MC is determined as a function of the photon energy. This correction is 
applied to the background component of qq7 MC events as well as to the hadronic side of the 
W~''W~7 MC events. The uncertainty on this correction is propagated to the measurement as 
a systematic uncertainty on the W"'"W~7 cross section which corresponds to ± 15 fb. 
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Determination of Anomalous Quartic Gauge Couplings 



The selected W"'"W~7 events allow us to constrain anomalous contributions to the SM quartic 
gauge boson vertices. In the framework of References El and the extended Lagrangian includes 
new dimension-6 operators, 

where oq/A^, ac/A^, and ctn/A^ are the AQGCs, and A represents the energy scale for new 
physics. The two parameters cto/A^ and ctc/A^, which are separately C and P conserving, 
generate anomalous W"'"W~77 and ZZ77 vertices. The term a„/A^, which is CP violating, 
gives rise to an anomalous contribution to W"'"W~Z7. Although there are already direct [7l|H] 
and indirect [12] limits on oq/A^ and Oc/A^, only the study of W+W~7 events allows for a 
direct measurement of the anomalous coupling a„/A^ through the W"'"W~Z7 vertex. 

The EEWWG program implements the effects of the AQGCs through the extended SM 
Lagrangian. Figure El shows how the anomalous coupling a„/A^ manifests itself through a 
deviation of the total cross section. 

The anomalous component from the above operators is linear in the photon energy at the 
matrix element level . This implies that also the shape of the photon spectrum is affected by 
AQGCs, in particular, the hard part of the energy distribution (see Figure Et)- The expected 
distribution for any value of the three AQGCs is obtained by reweighting each KORALW MC 
event with the ratio }V{E^, oq, a^, an) of the known differential distributions of at generator 
level: 

^^EEWWG j ^^KORALW 



yV(-E^,ao,ac,a„) = 77; (ao,a, 



The reweighting procedure is applied only to the ISR component of the MC selected sample, 
while the FSR (from KORALW) and the background components of accepted events are kept 
fixed. The possible dependence of the selection efficiency on the photon polar angle and on the 
angular separation from the charged fermions in the event is found to be negligible. 

Both the shape and the normalisation of the observed photon spectrum in the range from 5 
GeV to 30 GeV are used in a maximum-likelihood fit to each of the anomalous couplings ao/A^, 
a^/A^ and a„/A^, fixing the other two to zero. The effects of the same systematic uncertainties 
described for the cross section measurement are included, yielding the 68% CL intervals: 

-0.028 GeV"^ < ao/A^ < 0.028 GeV^^ 
-0.04GeV"2 < a^/A^ < 0.09 GeV'^ 

- 0.26 GeV"2 < a„/A2 < 0.23 GeV'^ 

The results are in good agreement with the SM value of zero for each of the anomalous quartic 
gauge couplings. The 1-parameter limits at 95% CL are: 

-0.045 GeV"2 < ao/A^ < 0.045 GeV"^ 
-0.08GeV'^ < a^/A^ < 0.13 GeV'^ 

- 0.41 GeV"^ < a„/A2 < 0.37 GeV'^ 



5 



The e+e — > 1/1/77 Process 

The selection of acoplanar multi-photon events is identical to that described in Reference EDI At 
least two photons with energies greater than 5 GeV and 1 GeV are required, with polar angles 
between 14° and 166°. The KORALZ [H] and NUNUGPV |22] Monte Carlo generators are 
used to model the e"'"e~ uu'j'j process according to the SM. The effects of the AQGCs Oq/A^ 
and ttc/A^ are simulated using the EENUNUGGANO program [Hj. Note that z/z/77 production 
is not sensitive to the fln/A^ coupling. 

We select 14 events at -y/i = 183 GeV and 21 events at i/i = 189 GeV compared to a SM 
expectation of 13.3 and 36.2 events respectively. 

The EENUNUGGANO program does not describe the effects of the SM s-channel Z ex- 
change diagrams and the interference between these diagrams and the W^W~ fusion diagram 
containing the W"'"W~77 vertex. Therefore additional cuts are applied to suppress the SM 
contribution. The energy of both photons must be greater than 10 GeV. If both photons are in 
the barrel region (| cos^| < 0.7), either the recoil mass must be less than 80 GeV or the sum of 
the photon energies must be greater than 100 GeV. If one or two photons are in the endcaps, 
where the SM contribution is larger, the recoil mass must be less than 75 GeV. After applying 
these cuts no data event is selected, consistent with the SM expectation of 0.15 events. 

The expected number of events for any AQGC value is calculated based on a sample of 
ten thousand simulated EENUNUGGANO events generated for several values of ao/A^ and 
Oc/A^. Its matrix element is used to reweight the events to any AQGC value required, testing 
the procedure by comparing the reweighted distributions to those from samples generated at 
various values of AQGCs. In all cases good agreement is observed. 

Since the program does not include higher order corrections due to ISR, these effects are 
estimated by implementing the EXCALIBUR |15j collinear radiator function. The cross section 
is reduced by about 16% which is used in the following. The remaining theoretical uncertainty 
of 5% ^3] is taken into account in the AQGC limits. The systematic uncertainty on the selection 
efficiency [20] gives a much smaller contribution. 

The 95% CL upper limit on the number of expected events from the AQGC signal is obtained 
taking into account the systematic error on the accepted cross section; this corresponds to the 
following 1-parameter limits at 95% CL: 

-0.067 GeV~2 < ao/A^ < 0.066 GeV"^ 
-0.18GeV-2 < aJA^ < 0.18 GeV'^ 

Conclusion 

All results obtained results are in a good agreement with the SM expectation of zero for each 
anomalous quartic gauge boson couplings. 

Combining the results on Oq/A^ and ad from our analyses of W^W~7 and z/i/77 produc- 
tion, we derive the following 1-parameter 95% CL limits: 

-0.043 GeV"^ < ao/A^ < 0.043 GeV'^ 
-0.08GeV"^ < aJA^ < 0.13 GeV'^ 
- 0.41 GeV"2 < an/K^ < 0.37 GeV'^ 
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Figure 1: Feynman diagrams containing a four-boson vertex leading to the (a) W+W 7 and 
to the (b) Ve^ell final states. 
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Figure 2: Differential distributions, for the semileptonic and fully hadronic W"'"W~7 decay 
modes, of (a) the photon energy, (b) the angle of the photon to the beam axis and (c) angle of 
the photon to the closest charged Icpton or jet. The hatched area is the background component 
from ZZ, Zee, and qq(7) events. Final state radiation includes the contribution of photons 
radiated off the charged fermions and photons originating from isolated meson decays. In the 
upper plot, the distribution corresponding to a non-zero value of the anomalous coupling a„/A^ 
is shown as a dotted line. 
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Figure 3: Measured cross section for e"*"e" W"'"W~7 at ^/s = 189 GeV (point) compared to 
the SM cross section as a function of the centre-of-mass energy (sohd hne) as predicted by the 
EEWWG Monte Carlo within the indicated phase-space cuts. The shaded band corresponds 
to the theoretical uncertainty of ±5%. The three dashed lines correspond to the cross section 
for non- vanishing values of the anomalous coupling a„/A^ (in GeV~^ units). 
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